Quantum-criticality and percolation in dimer-diluted 2D antiferromagnets 
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The S = 1/2 Heisenberg model is considered on bilayer and single-layer square lattices with 
couplings Ji, J2, and with each spin belonging to one J2-coupled dimer. A transition from a Neel 
to disordered ground state occurs at a critical value of gr = J2/J1. The systems are here studied 
at their dimer-dilution percolation points p* . The multi-critical point {g* ,p*) previously found for 
the bilayer is not reproduced for the single layer. Instead, there is line of critical points {g < g* ,p*) 
with continuously varying exponents. The uniform magnetic susceptibility diverges as T'" with 
a £ [1/2, 1]. This unusual behavior is attributed to an effective free-moment density ~ yi-"^ 'pj^g 
susceptibility of the bilayer is not divergent but exhibits remarkably robust quantum-critical scaling. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.40.Mg, 75.40.Cx 



A challenging and important aspect of quantum phase 
transitions [Ij is the influence of disorder (randomness) 
on the critical behavior, in the ground state as well as 
in the finite-temperature quantum-critical scaling regime 
0, For one dimensional quantum spin systems a 
real-space renormalization group (RG) scheme has 
rigorously established 0] a strong-disorder random sin- 
glet fixed point which controls the low-energy behav- 
ior for any strength of the disorder. This is an un- 
usual type of quantum criticality with dynamic exponent 
z = 00, which leads to, e.g., a susceptibility diverging as 
T^^ln~^(r). The RG procedure has also been carried 
out numerically for various two-dimensional (2D) mod- 
els, although in this case there is no rigorous proof of 
its validity. In the random transverse Ising model a ran- 
dom singlet phase was found 0], but only conventional 
critical points (finite z) were found for SU(2) symmetric 
(Heisenberg antiferromagnet) systems [7f. However, the 
disorder leads to unusual properties in the ordered and 
disordered phases 0, Q, e.g., quantum Griffiths effects. 

In the 2D = 1/2 Heisenberg model an order-disorder 
transition occurs at a critical strength of a coupling pat- 
tern favoring singlet formation on dimers, plaquettes, 
etc. [l, 0] . There is ample evidence from quantum Monte 
Carlo (QMC) studies [13, 111 that this transition is 
in the 3D 0(3) universality class, as predicted theoreti- 
cally 0, Il2| . According to the Harris criterion [l§], dis- 
order should be relevant at this transition. The applica- 
bility of efficient quantum Monte Carlo methods 14], in 
combination with a well developed RG scheme, makes 
dimerized Heisenberg models very well suited for ex- 
ploring disorder effects on quantum phase transitions. 
Recent studies have demonstrated a variety of scaling 
behaviors 0, [H, M, [13, El, [ll- The case of dilution 
disorder is particularly interesting, as it includes the spe- 
cial case of quantum-criticality at the classical percola- 
tion point, where geometrical and spin fluctuations are 
simultaneously divergent [l5, 16]. Here two cases of such 
multi-criticality are studied using QMC simulations — 
a bilayer and a dimerized single-layer model with ran- 



dom dilution of dimers, illustrated in Fig. [T] along with 
schematic phase diagrams. The dramatically different 
behaviors at the percolation threshold is the main focus 
of this Letter. The bilayer has a single point of quantum- 
criticality where the susceptibility scales to zero as T ^ 
In contrast, the single layer exhibits a line of 
quantum critical points where the susceptibility is diver- 
gent. 

The 2D Heisenberg model with random site dilution 
has been studied extensively because of its relevance to 
Zn doped cuprate antiferromagnets 2^, 21 1 . It was for a 
long time believed that site or bond dilution leads to a 
quantum phase transition before the geometric percola- 
tion point is reached 22, 2^ 24 1. However, recent QMC 
studies have shown that the transition coincides with 
the percolation point 25, |2^ and that the percolating 
cluster is ordered [2^. The critical exponents pertain- 
ing to equal-time correlations {v, (3, and ry) are therefore 
those of classical percolation. Other exponents {a,j,5) 
are given by combinations of percolation exponents and 
the dynamic exponent of the spin clusters j27| . 

The long-range order of the Heisenberg antiferromag- 
net on percolating clusters, which have fractal dimension- 
ality d = 91/48 [28(, implies that other couplings have to 
be introduced in order to realize a dilution-driven quan- 
tum phase transition [26| . This has been explored re- 
cently with the dimer-diluted bilayer model 15 , [13] • In 
the pure Heisenberg bilayer, with intra- and inter-layer 
couplings Ji and J2, there is an order-disorder transi- 
tion at g = J2/J1 ~ 2.52 [ll| (due to the tendency to 
singlet formation across the planes). If this system is 
diluted by randomly removing single spins, order is in- 
duced in the disordered phase because moments localize 
in the neighborhood of the vacancies. These moments 
interact and order at T = [2^. By instead diluting the 
two layers symmetrically, i.e., removing dimers consisting 
of nearest-neighbor spins on opposite planes, no local- 
ized moments form and a phase transition takes place at 
a coupling g which depends on the dilution fraction p, 
as shown in Fig. [TJ At the geometrical percolation point, 
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FIG. 1: (Color online) The lattices studied (left) and their 
T — phase diagrams (right). The couplings are JiSi • Sj 
and J2Si ■ Sj, with Ji, J2 > (in the bilayer all inter-layer 
couplings are J2 and in the single layer the vertical bonds 
alternate Ji,J2). Dimer dilution corresponds to removing J2- 
coupled pairs — such removed dimers are indicated by open 
circles. For the bilayer, the percolating cluster is ordered on 
the line {0 < g < g* jp"), whereas the single-layer is quantum- 
critical for {0 < g < g* ,p*). 



which clearly is the same as for a single site-diluted layer, 
p* « 0.41 (hole concentration) [28[, there is a multi- 
critical point {g* ,p*) at which the long-range order on the 
percolating cluster vanishes and the spins are quantum- 
critical. A critical coupling g* « 0.15 and dynamic ex- 
ponent z* ft! 1.3 were found in two independent QMC 
calculations 
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|16| |. The generic transition for p < p* 
has also been studied in detail by Monte Carlo simula- 
tions of an analogous 3D classical Heisenberg model with 
columnar defects 
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Also here an exponent 2; « 1.3 was 
found (but z* ^ z is expected because of the different 
cluster dimensionality). Long-range order in the pres- 
ence of quantum fluctuations on the line (g < g*,p*), 
which was believed not to be possible for a continuous 
order parameter [30| . was recently related to the fracton 
dimensionality of the percolating cluster [3l|] . 

The question now arises as to the generality of the be- 
havior found in the bilayer model. On its percolating 
cluster each spin has a neighbor in the opposite layer 
with which it correlates at low temperature. Magneti- 
zation fluctuations are thus quenched as T 0. Here 
the dimerized single-layer model is used to investigate 
the role of the bilayer symmetry upon dilution. Without 
disorder, the system has a quantum-phase transition, in 
the same universality class as the bilayer, at J2/J1 ~ 2.5. 
With dimer dilution again corresponding to random re- 
moval of J2 dimers, the percolation point p* — 1/2 be- 
cause the dimers are connected as a triangular lattice [2^ . 
The T = phase diagram, outlined in Fig.[Tl is similar to 
that of the bilayer in that there is a flnite segment of the 
phase boundary &i p = p* , terminating at a point {g* ,p*) 
beyond which the transition occurs for p < p* . However, 
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FIG. 2: (Color online) Temperature scaling of the suscepti- 
bility of the infinite percolating bilayer cluster for different 
coupling ratios g, using a dynamic exponent z — 1.36. 



as will be discussed in detail below, there is a striking 
difference: Whereas the percolating cluster of the bilayer 



has Neel order for < g < g* [15|, \lu\ , the percolating 
cluster of the single dimerized layer is quantum-critical on 
the whole line (0 < 5 < g* ,p*), with g* « 1.25. On this 
line, the magnetization fluctuations are not completely 
quenched as T — > 0, leading to a divergent susceptibil- 
ity; X ^ T^°', with a — + 1/2 for g g*^ and a ^ 1 
for g —^ 0. For g = 0, a, Curie susceptibility is expected 
on account of the percolating cluster breaking up into 
smaller pieces when all the J2 couplings vanish (the per- 
colation point here is Pq » 0.29). Some of these clusters 
contain an odd number of spins. For g > 0, the divergent 
susceptibility can then be attributed to effectively iso- 
lated subclusters with net moments, which are gradually 
"frozen out" as the temperature is reduced. The form 
X ^ corresponds to a free-moment density scaling 
as T^^" . This remarkable behavior will here be demon- 
strated on the basis of large scale QMC (stochastic series 
expansion [lj|) calculations. Only results exactly at the 
percolation threshold, p = p*, will be discussed. The 
numerical techniques and special methods developed for 
studies of random systems at ultra-low temperature are 
discussed in detail in [2^ . 

The temperature dependence of the uniform suscepti- 
bility x{T) of the bilayer close to the multi-critical point 
was discussed before in Ref. flBj (averages over all clus- 
ters were presented in [IGi]). Fig.[2]shows a more extensive 
set of high-precision results for the largest cluster on LxL 
lattices with L = 256 {L —^ 00 converged for the temper- 
atures shown) at p* = 0.4072538 [S^]. Averages over sev- 
eral thousand dilution realizations were taken. The tem- 
perature is scaled according to the expected quantum- 
critical form, X = 0, + foyi/^-i Q ^ where a and b are con- 
stants and a = at a quantum-critical point. Using the 
fractal dimension d = 91/48 and adjusting z to obtain 
a linear x versus T'*/^-!, the same dynamic exponent. 
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FIG. 3: (Color online) Temperature dependence of the suscep- 
tibility of the percolating single-layer cluster, (a) At and close 
to the coupling g* « 1.247. (b) T"^''^ scaling for 1 < gr < p*. 
(c) Scaling with varying exponent for g <\. 



z = 1.36±0.01, is found for all < g ^ g* . An improved 
estimate g* = 0.118 ± 0.006 is also obtained. Note that 
bilayer criticality can be expected only for T < g, which 
is indeed the case in Fig. [2] 

In a clean quantum-critical system, there is a low- 
temperature cross-over of x(T) to a "renormalized classi- 
cal" behavior when g < gc, at a temperature of the order 
of the spin stiffness ps [3]. No such cross-over is seen 
in Fig. [21 however. Although it cannot be completely 
excluded that there is a cross over at still lower temper- 
ature, one can also argue that there should be no cross- 
over, because = on the percolating cluster 1^ 2^ 
(although there is long-range order — this unusual behav- 
ior has also been discussed in [s^l ) and there is no appar- 



ent energy scale, except T, to govern the long-distance 
physics of the spins on this fractal network. Thus the 
multi-critical point {p*,g*) may control the T < g tem- 
perature dependence for all < g < g* at p* . Consider- 
ing that the scaling in Fig. [2] extends down as low as to 
T = J/256 and to couplings g ^ g* /3, the results appear 
to support such an unusual manifestation of quantum- 
criticality on a fractal percolating cluster. 

Turning now to the single-layer model at its dimer per- 
colation point, p* = 1/2 for g > 0, the disorder-averaged 



uniform susceptibility of the largest cluster was calcu- 
lated for L up to 256, down to T = J/512 {L ^ oo 
converged). Here a special point {g*,p*) is found which 
separates qualitatively different behaviors of the uniform 
susceptibility. As seen in Fig.EJa), at g* = 1.247±0.001, 
X is linear in T and approaches a finite value as T — > 0. 
For g > g* the susceptibility drops to zero and for g < g* 
it diverges. As shown in Figs. mb),(c), the divergence 
is of the form x ~ T^", with a very close to 1/2 for 
1 ~ < 5* a —^ 1 for g —^ 0. Such a divergence 
can be interpreted as a temperature dependent fraction 
^ rpi-a q£ effectively free magnetic moments. 

As already noted, exactly at 5 = the percolating 
cluster is broken up into smaller subclusters and then 
a Curie behavior, x ^ T^^, is expected on account of 
clusters with a net spin. For small but non-zero g, one 
might then have expected a cross-over from Curie be- 
havior when T > g to a finite susceptibility as T — *■ 0. 
Instead, it appears that coupling the subclusters leads 
to a g dependent power-law temperature scaling of the 
number of effectively free moments. The effective cou- 
plings of these moments to each other must thus have 
a g dependent power-law distribution, leading to a self- 
similar structure of free moments different from that of 
the underlying fractal cluster. It is remarkable that this 
behavior persists also when 5 « 1, where the picture of 
weakly connected subclusters is not obviously relevant. 

The bila yer percolating cluster is ordered at T = for 
9 < 9* Si- This is not the case for the single dimer- 
ized layer at p* . Instead, quantum-critical fluctuations 
are observed for Q < g < g* ■ Consider the staggered 
structure factor S{ti, tt) and susceptibility xC""? of a 
cluster of Nc spins, 

1 

K 



Y(7r, tt) — — 



dr{Stir)SUO)), 



where = 



(^^I'jXi+vi-xj-yj ^ At a quantum-critical 
point, these quantities, averaged over disorder, should 
scale with the system size as {S{TT,n)/Nc) ~ 
(x(7r, 7r)/iVc) ~ L'O', with 7^ = -{z + 77) and 7^ = -ry 
(normalizing by Nc ~ L'' before disorder-averaging leads 
to some reduction of statistical fluctuations). Fig. [3] 
shows results for g = 1. The observed scaling gives 
z w 3.0 and rj « —1.7. The dynamic exponent can be 
compared with the expected quantum-critical suscepti- 
bility; X T-" with a = 1-d/z^. With the exponent 
a = 1/2 obtained from x(T) in Fig. [3]Jb), it is apparent 
that this relationship does not hold here [z extracted from 
S{7T, tt) and x(7r, tt) should be the actual dynamic expo- 
nent]. At the special point g* the extracted a = and 
z Ki d (not shown here) are in fact consistent with this 
relationship. This is also the case at the bilayer multi- 
critical point 15|. The single-layer quantum-criticality 
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FIG. 4: (Color online) Finite-size scaling of the staggered 
structure factor and susceptibility, normalized by the cluster 
size, at intra-dimer coupling g = 1. The lines show scaling 
with exponents indicated in the legends. The results where 
obtained at T sufficiently low to give the ground state. 



on the line (0 < g < g*,p*) thus appears to be funda- 
mentally different. It should be noted that the behavior 
is not consistent with a Griffiths phase |7], as the spin 
correlations in that case should be exponentiaUy decay- 
ing, in contrast to the power-law seen in Fig. 2) 

An extended line of quantum-critical points was not 
anticipated on the basis of a real space RG approach 
developed recently for quantum rotors on a percolating 
cluster 



3l|. The critical line discovered here for the 



single-layer model is more similar to the ID random sin- 
glet phase [51], in that there is a temperature dependent 
fraction of effectively free moments. In the random sin- 
glet phase z = 00 whereas in the model studied here z 
is finite and diverges in the limit g 0. A temperature 
dependent fraction of effective moments has also been 
observed in a 2D model of interacting localized moments 
[3^ . However, there the asymptotic T — > susceptibil- 
ity is always Curie-like, and there is long-range order at 
T — 0. An RG calculation for random frustrated mo- 
ments in the continuum shows a x{T) divergence with 
varying exponent [s^ . The ground state properties were 
not accessible in that study. 

The bilayer multi-critical point {p*,g*) has been ar- 
gued [1^ [3| to influence finite-temperature properties 
of single-layer Zn doped cuprate antiferromagnets, for 
which a dynamic exponent z w 1.4 was found in neu- 

However, although g* 



tron scattering experiments 
is small (« 0.12) it is difficult to explain how bilayer 
quantum-criticality could be realized when T ^ g = 
(due to an expected cross-over at T « 5; see Fig. [2]). 
Physical realizations of single-layer dimer dilution are 
not immediately obvious. Nevertheless, the results pre- 
sented here serve to illustrate rich and surprising behav- 
iors arising from the interplay of classical percolation and 
quantum fluctuations, going beyond previous examples of 
scaling in percolating fractal structures (s^l . 

This work was supported by the NSF under Grant 
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